Echolocation by Quasiparticles 
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It is shown that the local density of states (LDOS), measured in an Scanning Tunneling Microscopy (STM) 
experiment, at a single tip position contains oscillations as a function of energy, due to quasiparticle interference, 
which is related to the positions of nearby scatterers. We propose a method of STM data analysis based on this 
idea, which can be used to locate the scatterers. In the case of a superconductor, the method can potentially 
distinguish the nature of the scattering by a particular impurity. 
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Scanning Tunneling Microscopy (STM), which measures 
the "local density of states" (LDOS) as a function of position 
and energy set by the bias voltage, has opened the door to 
imaging the sub-nanoscale topography and electronic struc- 
ture of materials, including normal metals and especially 
cuprate superconductors |2ll3l|4l|J,|aLZl|8|,|9|]. 

The dispersion relations of (Landau or Bogoliubov) quasi- 
particles may be extracted from STM data on normal met- 
als ifToifTlll and superconductors lll3ll . via the inverse method 
called Fourier transform scanning tunneling spectroscopy 
(FT-STS) CSEl, or directly in real space Ol. This tech- 
nique is based on the fact that impurities produce spatial mod- 
ulations of the LDOS in their vicinity - standing waves in 
the electronic structure that generalize the Friedel oscillations 
found in metals at the Fermi energy. In the cuprates BSCCO 
and CaCuNaOCl lfl3ll . experiments showed these quasipar- 
ticle oscillations were dominated by eight wavevectors that 
connect the tips of "banana" shaped energy contours in recip- 
rocal space, the so-called Octet model as explained theoret- 
ically [12]. For optimally doped samples, the dispersion in- 
ferred from these wavevectors agrees well with d-wave BCS 
theory indicating the existence of well-defined BCS quasipar- 
ticles in this regime. 

The central observation of this paper is that the same 
Friedel-like oscillations of the LDOS, analyzed in the 
space/momentum domain by FT-STS, are also manifested in 
the energy/time domain. Our analysis shows that the small 
impurity-dependent modulations of the LDOS have a period, 
in energy, inversely proportional to the time required by a 
quasiparticle wavepacket to travel to the nearby impurities and 
back - hence we call it "quasiparticle echo". From this, in 
principle, one can determine the location and (in a supercon- 
ductor) the nature of the point scatterers in a particular sample. 

Quasiparticle echo — The basic idea of the LDOS mod- 
ulations may be understood semiclassically. The LDOS 
N(r;u>) is defined as — (l/7r)ImG(f, r; u>), the time Fourier 
transform of the local (retarded) Green's function G(r, r; t). 
Imagine a bare electron wavepacket (centered on energy uS) is 
injected at time t = at point r in a two-dimensional mate- 
rial: the Green's function expresses its subsequent evolution. 
Assuming there are well-defined quasiparticles at this energy 
with dispersion E(k); then for every wavevector k on the en- 
ergy contour E(k) — u>, the wavepacket has a component 



spreading outwards at the group velocity v g (k) = Vg-E'(fe) / H. 
When this ring reaches an impurity at r*i mp , it serves as a 
secondary source and the reflected wavepacket arrives at the 
"echo time" 
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for the k such that v g (k) || R = r; mp — r. This creates a sharp 
peak at t = T e in G(r, r; t) [see Fig.[T](d)], and hence modu- 
lations as a function of uj in its Fourier transform N(r; ui) with 
period Auj = 2irh/T e Generically, for a particular im- 
purity direction, \v g | varies with energy, so the the modulation 
in SN(lu) due to the impurity is "chirped" correspondingly. 

We illustrate the quasiparticle echo first by a numerical cal- 
culation for a normal metal, defined by the lattice Schrodinger 
equation for the wavefunction Ui on site i: 
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(Uj + = Eui 



(2) 



Here the t'^s are intersite hoppings and the /i,'s are on-site 
potentials (including the chemical potential); in this paper, 
we assume they are translationally invariant except at dis- 
crete (and dilute) impurity sites. We take the specific case 
of nearest-neighbor hopping t at half-filling, so the the dis- 
persion is e(k x ,k y ) = ~2t(cosk x + cosfcj,), and we place 
one (repulsive site potential) impurity at the origin. To nu- 
merically calculate the LDOS, we used the Recursion method 
lfl6ll . which is well-suited for cases without translational sym- 
metry. 

Fig- E a ) shows the impurity case LDOS which has echo 
oscillations on top of what otherwise would have been clean 
case LDOS, visible along the sides of the peak. Note that, 
for us to see more than one oscillation within the bandwidth, 
the impurity must be at least several sites away; hence the 
oscillations always have small amplitude and are best viewed 
by subtracting the clean LDOS. Throughout the paper, energy 
is in units of t and time in units of t^ 1 with t = 1 and H = 1. 

For a given energy ui, we define Aw(w)/2 as the separa- 
tion of the zeroes that bracket u> in the (subtracted) SN(lu) 
trace, and let T e (uo) = 2%H/ Aw(w), We chose E = 0.7t 
and R in the [1,1] direction, for which the group velocity is 
v g = 2.785t/h. Then, using 6N(20, 20; w), <W(30, 30; uj), 
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FIG. 1: LDOS as a function of energy, showing oscillations due to 
quasiparticle echoes, (a) LDOS at a point 20\/2 away from a point 
impurity along the [1, 1] direction(lattice constant = 1). (b) corre- 
sponding LDOS after subtracting the clean LDOS : 8N(20, 20; to), 
(c) <5JV(40, 40; uj), (d) Magnitude of Local Green's function as a 
function of time : jG(20, 20; T)\. The singularity appears at time 
T e /2, where T e is given by iff). As we change the distance along this 
direction, the shortest echotime changes in proportion in accordance 
with our semiclassical expectations. 




FIG. 2: (a) Schematic of few measurements around an impurity. The 
arrowheads represent the STM Tip postions. After measurement, we 
get (b) 5N(f A ;uj), (c) 8N(f B ;uj) and (d) SN(f c ;uj). Extracting 
the echotimes for each measurement at uj — 0.7t, we locate the im- 
purity, shown as a black dot, in the first panel. Note, that the locus of 
impurity locations changes with uj, and is of the shape shown only at 
uj = 0.7t. 



and 6N(40, 40; uS) [the first and last trace of these are shown 
mFig.Ecd)], we read off Aw/2 = 0.1545, 0.103, and 0.077, 
from which v g T e /2 = 20.04^, 30.05^2, and 40.22a/2, re- 
spectively. The proportionality between the oscillation rate 
and the actual distance confirms the semiclassical explanation 
of these modulations. 

Echolocation — Using these quasiparticle echoes, we can 
locate the position of impurities by measuring the LDOS wig- 
gles at a few points in the vicinity. At each point, we ex- 
tract the wiggle period Auj and hence the echo time T e = 
2tt/Alu. Then (Q]) defines a locus of possible impurity lo- 
cations, {v group (k)T e /2 : e(k) = uj. The intersection of 
the loci from STM spectra taken at multiple points r will lo- 
cate r; mp uniquely. Furthermore, via a more exact derivation 
of the LDOS modulations (see below), the amplitude of the 
LDOS modulations tells the scattering strength of the impu- 
rities (in Born approximation they are proportional to each 
other). Once an impurity has been pin-pointed, the higher- 
energy STM spectrum at that point may independently iden- 
tify the chemical nature of the impurity, e.g. in cuprates 11511 
and thus may reveal which kinds of impurities are important 
for the scattering of quasiparticles. 

As a test, we evaluated the subtracted LDOS at three points 
r A = (-30,0), f B = (-20,20), and r c = (15,30), with 
the impurity at r = 0. From the half-periods of the wiggles at 
energy = 0.7i, (extracted as before) we found the respective 
echo times T A = 39.9, T B = 20 .4 and T c = 36.7. The three 
scaled loci(scaled by half the respective echotimes), shown in 
Fig. |2]e), intersect at (0,0) as can be seen graphically, thereby 
demonstrating the idea of echolocation. A more careful nu- 
merical analysis can be done to extract errors in echolocation 
as well. 



Analytic derivation — Adopting the T-matrix formalism, 
we can obtain an analytic form for the LDOS modulations. 
Formally, the difference in dirty LDOS and clean LDOS for a 
single point impurity is given by 
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G (r- f imp ; w)T(w)G (r il 
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where G (r,f imp ;w) = G {r - r imp ;w) = G (R;lo) is the 
free propagator; LDOS modulations are due to interference 
between the two Go factors. 



Gq(R;lo) = lim 
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^° + Jb.z. (2tt) 2 cj + iS-e(k) 
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The integrand is singular all along the energy contour e(k) = 
uj, which we also parametrize as k e (s), where s is the arc- 
length in reciprocal space. By the change of variables z = 
e %k v we conver t the inner (k y ) integral to a complex contour 
integral in the z plane (rewriting e(k x , ky) as an analytic func- 
tion of z); for k x values found on the energy contour, the z 
path encounters two poles, one inside and one outside, de- 
pending on the sign of S. Extracting the residue and absorbing 
factors, we get 
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where 77 (s) = 1 on the half of the energy contour where 
sgn(<5) = sgn(\v g (oj, s)\) and zero on the other half. The non- 
singular term G n0 n-singuiar comes from the integrals over k y 
which do not cross the energy contour. 

At large R, the two-dimensional BZ integration will be 
dominated by those k 111 811 on the energy contour where the 
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phase in the numerator is stationary, i.e. v g (k) || R: let us call 
such a point kg (so it is a function of the direction R and of 
lo). Using standard formulas of the stationary phase approxi- 
mation IU9I1 we get asymptotically 
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Here k 1 is the curvature d?k e /ds 2 of the energy contour at 
Using (Q~|i and (0, we finally get 
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cos(2k & (R,u))-Ry (7) 



valid in the limit of a distant impurity. (All factors are actually 
functions of R and lo: these arguments are shown only in the 
rapidly varying factors.) As we change lo to lo + Slo keeping R 
fixed, the chain rule gives kg(io+Suj) — kg(R, lo) = v 



- l 8ioR 



so, with <f> = ■ R, we get 



cos hk A (R, lo + 5lo) ■ Rj -> cos(0 + TJlo). (8) 

This confirms the simple semiclassical prediction Alo = 
2n/T e (see Eq. ([TJ) for the modulation period due to echoes. 
The same quasiparticle interference is responsible for the spa- 
tial oscillations evident in (O and the energy oscillations in 
©. 

Echoes in cuprate superconductors — Additional relevant 
issues arise in case of superconductors. To discuss these, we 
use a mean-field Bogoliubov-DeGennes(BDG) Hamiltonian 
with/without a single point impurity as shown below. 
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where we are using a lattice formulation of BDG equations. 
The u.i& and w^s represent particle and hole amplitudes on site 
i, UjS and piS represent the intersite hoppings and site chemi- 
cal potentials respectively, and Ay represent the off-diagonal 
order parameter amplitude. We discuss d-wave superconduc- 
tors (dSC's) to highlight this method's application to cuprates. 
For dSCs, Ay is nonzero only on nearest-neighbor bonds and 
A- - , - = — A- - , - because of the d-wave nature. Our normal 
state is the same nearest neighbor tight binding model on the 
square lattice with t = 1 and off-diagonal hopping amplitudes 
set to |A| = 0.1. The Recursion method was extended to su- 
perconductors in ifnll and is used for our numerics. In Fig. 
[3]c) and d), we show the LDOS(after subtracting the clean 
LDOS shown in Fig. 0a)) at 20^2 distance from an impurity 
along the (1,1) direction for the case of a potential scatterer 
and an anomalous pair potential scatterer (which scatters an 
electron into a hole and vice versa) respectively. 

In contrast to the normal case, there are two different wig- 
gles : a fast one and a slow one. The reason for this is that the 
dSC quasiparticle dispersion gives rise to two different group 
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FIG. 3: Quasiparticle echoes in a d-wave SC. (a) no impurity 
No(lo) showing the d-wave gap, (b) A caricature of the two dif- 
ferent group velocities along (1,1) direction for d-wave Bogoliubov 
dispersion, (c) <5iVoRD(20, 20; lo) for an ordinary impurity and (d) 
5JVano(20, 20; lo) for an anomalous impurity. 



velocities in the (1,1) direction M20I1 . We also note that the 
fast wiggles exist only within the gap while the slow wiggles 
are both inside and outside the gap. In Fig. [3]b), we show 
the constant energy contours for the quasiparticle dispersion 

given by E{k) — \J e(fc) 2 + A(fc) 2 , the gradient of which is 
the quasiparticle group velocity. From Fig. [3]b), we see that 
along (1,1), the banana-shaped energy contours in the first 
and third quadrants give one velocity (which corresponds to 
the slow wiggles), while the contours in the second and fourth 
quadrants give a slower velocity (which corresponds to the 
fast wiggles). For E > | A|, there are no longer "banana" con- 
tours, so we get only one group velocity (similar to the normal 
case) and hence only one kind of wiggle is seen in Fig.0c,d) 
outside the cusps. 

Once the impurity is located using the loci intersection 
method desribed before, one can study the LDOS data around 
the impurity to infer the impurity's strength and whether it is 
ordinary (magnetic/nonmagnetic) (cf . Ref . |2l| and references 
therein) or anomalous M22I1 . This distinction is already visible 
in individual spectra: provided the normal state is particle- 
hole symmetric, one gets particle-hole symmetric echo oscil- 
lations SNaho from an anomalous impurity, since it scatters 
electrons into holes and vice versa [Fig. [3d)]; this is not the 
case for (Word from an ordinary impurity [Fig.fSfc)]. 

A second diagnostic distingushing (nonmagnetic) ordinary 
scatterers from anomalous ones is the real-space pattern of the 
surrounding standing waves in the LDOS, which is best seen 
in Born Approximation. In this limit, the impurity T-matrix is 
of the form (in the 2x2 Nambu notation) U- lmp Ts or Ai mp Ti 
for the ordinary or anomalous cases, respectively. Then the 
echo oscillations take the respective forms 



5N 0RD oc U imp (G 2 n -G 2 u 
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id (2G11G12). 
(10) 
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Here, the G^ s are the matrix elements of the usual free prop- 
agator G (fc;w) = - ^(fc) 2 )- 1 ^ + e(fc)r 3 + A(fc)n] 
thus in real space 

G (R;u) = J f )( s )y5(M s ! w )) + Gnon-singuiai- (11) 

where g(p; R,l>) = 1+ ~(e(p)T 3 + A(p)n). 

We can carry out the stationary phase approximation as be- 
fore, but instead we numerically calculated the propagator us- 
ing Eq. ( TTTT i. since we are interested in LDOS information 
around(close) to the impurity. In Fig. |4] we show SN around 
an impurity over a grid of 20x20 lattice points(shown one 
quadrant with others related by symmetry). 




FIG. 4: Shown are the SN(R; lo = 0.35t) around an impurity over 
a grid of (0,20)x(0,20) with other quadrants related by symmetry, 
(a) <5-/Vord, (b) SNano- A subgap value of ui — 0.35t was chosen 
arbitrarily. 

We see that certain of the real-space oscillations, present in 
case of the ordinary impurity, are suppressed in the case of 
a d-wave anomalous impurity. This is the same effect as the 
suppression of certain "octet" vectors for the case of 

d-wave anomalous impurity as argued in [22]'s Eq. 10 and the 
following paragraph. Our real-space analysis qualitatively du- 
plicates that of Refl22l illustrating how the real-space QPI and 
our energy-domain echoes are complementary manifestations 
of the same phenomenon. 

Conclusion and Discussion — In conclusion, we have in- 
troduced a method of STM data analysis in the energy do- 
main as a phenomenological tool for the study of real materi- 
als, complementary to FT-STS. Since it is based on the same 
quasiparticle interference effects already used successfully in 
FT-STS, we have confidence that the signals will be observ- 
able. They should be particularly strong in materials with an 
energy-dependent group velocity in some range of energies, 
such as d-wave superconductors and also graphene II23I1 . 

Since the echo analysis can be done in local patches of the 
sample (unlike FT-STS which fourier transforms over a larger 
region), we can locally verify the existence of quasiparticles 
at various energies through QPI. In particular, in cuprates, 
echoes might be used to check the hypothesis of quasiparti- 
cle extinction 12411 above a certain energy. Furthermore, we 
have argued that echo analysis might reveal the nature of spe- 
cific impurities 12511 in a sample, information which hitherto 



was (at best) known statistically. 
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